Theory of decoherence in a matter wave Talbot-Lau interferometer 
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We present a theoretical framework to describe the effects of decoherence on matter waves in 
Talbot-Lau interferometry. Using a Wigner description of the stationary beam the loss of interference 
contrast can be calculated in closed form. The formulation includes both the decohering coupling to 
the environment and the coherent interaction with the grating walls. It facilitates the quantitative 
distinction of genuine quantum interference from the expectations of classical mechanics. We provide 
realistic microscopic descriptions of the experimentally relevant interactions in terms of the bulk 
properties of the particles and show that the treatment is equivalent to solving the corresponding 
master equation in paraxial approximation. 
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I. INTRODUCTION 



The art of demonstrating the wave nature of mate- 
rial particles experienced considerable advances in recent 
years; see 0, 0, Q an d references therein. The interfer- 
ing species evolved from the elementary particles of the 
early experiments |J. [of to composite objects with an in- 
ternal structure. In particular, the experiments in atom 
interferometry have left the stage of proof-of-principle 
demonstrations, and provide substantial applications in 
metrology 00, 00- Objects with even larger complex- 
ity, such as molecules or clusters, exhibit a rich internal 
structure that can interact in various ways with exter- 
nal fields. Their interference is highly sensitive to the 
corresponding phase shifts, thus offering the potential to 
measure molecular properties with unprecedented preci- 
sion. At the same time any coupling to uncontrollable 
fields and environmental degrees of freedom severely lim- 
its the ability of large objects to show interference. These 
effects are bound to become dominant as the chosen ob- 
jects increase in size and complexity. 

The influence of environmental coupling on a quantum 
system may be described by decoherence theory ^) . 
It considers both the influence of noise due to uncontrol- 
lable external fields and the effect of the entanglement 
with unobserved dynamic degrees of freedom. This lat- 
ter phenomenon - the dynamic derealization of quantum 
coherence into many environmental degrees of freedom - 
largely explains the emergence of classical behavior in 
a quantum description. In particular, it describes the 
wave-particle complementarity encountered if one seeks 
to determine by a (macroscopic) measurement device the 
path taken. Since matter wave interferometers establish 
quantum coherence on a macroscopic scale they are sen- 
sitive tools to probe the quantum-to-classical transition 
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of complex objects. 

The purpose of this article is to provide the theoret- 
ical framework needed to describe the diffraction and 
decoherence effects encountered in the interferometry of 
large, massive objects. We focus on near-field Talbot- 
Lau interference, which is the favored setup for short de 
Broglie wavelengths. We take care to describe the effects 
of diffraction and decoherence with realistic parameters, 
to permit a direct quantitative comparison with the ex- 
perimental signal. The interactions are treated on a mi- 
croscopic level using the bulk properties of materials and 
particles. We note that the recent interference experi- 
ments with fullerenes and biomolecules 0, 0, 
were analyzed using the theory presented in this article. 

Before going into calculations we start with an informal 
discussion of Talbot-Lau interference 0, 0, 0] . In this 
setup an essentially uncollimated particle beam passes 
three parallel gratings. Effectively, the first grating acts 
as an array of collimation slits which illuminate the sec- 
ond grating. Diffraction at the second grating then leads, 
for particular choices of the grating periods and the wave- 
length, to a high contrast near-field interference pattern 
at the position of the third grating. This density pattern 
is observed with the help of the third grating by recording 
the transmitted flux as a function of the lateral grating 
position. 

An important advantage of the Talbot-Lau effect is the 
favorable scaling behavior with respect to larger masses 
of the interfering object |l9|. Unlike in far-field diffrac- 
tion, where the required grating period falls linearly with 
the de Broglie wavelength, it decreases merely like the 
square root in the Talbot-Lau setup. In addition, the 
collimation requirements are much weaker than for far 
field diffraction, and the spatially resolving detector is 
already built into the device. 

However, for a fixed particle velocity it is not imme- 
diately evident whether the observed signal proves gen- 
uine quantum interference, since a certain fringe pattern 
could also be expected from a classical moire effect. This 
classical pattern can be suppressed by an appropriate 
choice of the open fraction of the grating and, unlike the 
strong wave-length dependence found in the interference 
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effect, the ideal classical shadow fringes do not depend on 
the velocity. Nonetheless, in order to distinguish clearly 
the quantum phenomenon from a classical expectation it 
is necessary to describe the quantum and the classical 
evolution in the same theoretical framework, thus ensur- 
ing that all interactions and approximations are treated 
equally. 

A first aim of this article is to provide such a descrip- 
tion that draws the demarcation line between the predic- 
tions of quantum and classical mechanics concisely and 
quantitatively. The second aim is then to account for 
the relevant environmental interactions, thus providing a 
quantitative description of the transition from the quan- 
tum to the classical behavior. For both goals it will be 
helpful to describe the state in the interferometer in terms 
of a stationary, unnormalized Wigner function. Due to 
the stationary formulation the effect of decoherence will 
not be given by a master equation. Thcrfore, it is shown 
in the final part of the paper that our treatment is equiv- 
alent to the conventional dynamic formulation of deco- 
herence in terms of a normalized Wigner function. 

The structure of the article is as follows: In Sect.lTTlwe 
review the coherent Talbot-Lau effect and give a formula- 
tion in terms of the Wigner function. The corresponding 
classical shadow effect is calculated on an equal footing 
in the phase space representation. The influence of the 
interaction with realistic gratings, which is very impor- 
tant for a quantitative description, is accounted for in 
Sect. 11111 Those effects are also treated on an equal de- 
gree of approximation in the quantum and the classical 
description. In Sect. II VI we include the possibility of de- 
coherence and show how it can be accounted for analyt- 
ically. The specific predictions for decoherence due to 
collisions and due to heat radiation are then obtained in 
Sect. [V] In Sect. IVII we relate the description of deco- 
herence in terms of a stationary beam to the solution 



of the corresponding time-dependent master equation. 
Concluding remarks are given in Sect. IVlH 



II. THE TALBOT-LAU EFFECT IN THE 
WIGNER REPRESENTATION 

Since the coherent theory of the Talbot-Lau effect can 
be found in the literature [2(1 l2ll l22l l23| we shall present 
no detailed derivations, but discuss the approximations 
involved and state the results in terms of the Wigner 
function as far as needed for the later inclusion of deco- 
herence effects. Consider the usual interferometric sit- 
uation where a flux of particles enters at z = with a 
longitudinal momentum p z that is much greater than its 
transverse components. Ideally, the particle is in a mo- 
mentum eigenstate, or in an incoherent mixture thereof, 
before passing a number of collimation slits and gratings. 
The vector r = (x, y) describes the distance of the par- 
ticles from the interferometer axis. In the usual paraxial 
approximation this separation, as well as the structures 
in the grating and in the collimation planes, are assumed 
to be small compared to the distances Li between the 
optical elements, |r| < L,. In this case one may evaluate 
the transmission to leading order in |r|/Z,-. This approx- 
imation implies that the longitudinal and the transverse 
part of the state remain separable throughout the inter- 
ferometer. It follows that the discussion may be confined 
to the transverse degrees of freedom as described by ?/>(r) 
if the evolution is completely coherent — or, in the gen- 
eral case, by the density matrix p(r,r'). 

Given the wave function V'o( r ) on the z = plane the 
free unitary evolution up to the plane z = L yields, to 
leading order, the transverse state, 



^oex P (i!^)*oM +£>(£) , (1) 

as follows from an asymptotic expansion of the free Green function, e.g. [24|. An important feature of this paraxial 
approximation is the fact that it reflects the composition property of the exact propagation without any loss of 
accuracy. That is, propagating the wave function first by a distance L\ and subsequently by the distance L2 yields 
exactly the same result as a single propagation by L\ + L2, 



2nta(Li+L 2 ) J \ h2(L 1 + L 2 

r 



which follows from Gaussian integration. Hence, within tervals and integrating over the interjacent planes. This 
the paraxial approximation no loss of accuracy is intro- freedom will be used below as a crucial ingredient when 
duced by dividing the propagation into a sequence of in- we describe the effects of decoherence. 
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Note that the composition property © does not re- 
quire a large separation between the planes. Even for 
infinitesimally close planes one obtains the correct ex- 
pression QJ. This can be seen immediately in J2J by 
noting a particular representation of the two-dimensional 
S function jH eq. (A.33)], 



Pz 



2-KhiL 



exp 1 



■ Pz 



2L 



^(r-r ) 



(3) 



The first equality in J2J) also shows how the existence 
of an ideal grating at z = L\ would affect the propagated 
state. In the case of a binary grating the reintegration 
would be simply restricted to the transparent parts of 
the grating plane. In general, an ideal grating causes an 
amplitude and phase modulation 

^l(r 2 ) =i(r 2 Wi(r a ) with |t(r 2 )| < 1 , (4) 

which is accounted for by the additional appearance of a 
grating function i(r 2 ) under the integral. 

The passage of a particle stream through a general 
interferometer may now be described as a sequence of 
transmissions through gratings or collimation slits as de- 
scribed by 10} each followed by a free evolution . This 
holds also for general, mixed states since any density op- 
erator can be represented as a convex sum of projectors 
to pure states. 



A. Wigner function 

We proceed to formulate the propagation in the 
Wigner representation, which has several advantages. 
First, the Wigner function permits a direct comparison 
of the quantum evolution to the classical dynamics in 
terms of phase space distributions. Second, and more 
importantly, it is the most convenient starting point to 
include effects of decoherence in Sect. IIVI Finally, the 
free evolution has a particularly simple form in the 
Wigner representation. 

The Wigner function is the Fourier transformation of 
the position density matrix p(r, r') with respect to the 
two-point separation A = r — r' 26] , 



w(r, p) 



1 



dAe i P AA 



A 

y 



p \ v ?r i r i 77 



(5) 

It may be viewed as a quantum analogue to the classical 
phase space distribution /(r,p), with p the transverse 
momentum vector. 

In order to obtain the free unitary evolution of the 
Wigner function we note that the density matrix in po- 
sition representation has the general form 



p(r,r') = J d M5 ( M )^(r)^(rO 



(6) 



with J d/ig(/i) = 1. According to P a free unitary evo- 
lution by the distance L yields 



P(r,r') = 



(2irh) z L 2 



drodi-Q exp ( i 



| r-r | 2 - l^-r^l 2 
2L 



(7) 



From JSJ and J7J it follows that a free unitary evolution 
by the distance L changes the Wigner function according 
to 



WL(r,p) = w r 



— P,P 

Pz 



(8) 



This transformation is particularly simple and, as one 
expects, is identical to the free movement of a classi- 
cal probability density in the phase space of the trans- 
verse degree of freedom. The decisive difference between 
the classical and the quantum phase space dynamics is 
found in the transformation for passing through a grat- 
ing. Equation Q implies that by going through a grating 
the Wigner function undergoes a convolution 



w'(r,p)=J dqT(r,q)w(r,p - q) 



(9) 



which in general builds up quantum coherences that show 
up as oscillations in the momentum direction. Here we 
define the convolution kernel analogously to (0) as 



T(r,p) 



(2nhy 



dA e ipA/ri i 
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~2 



t* 



A 

~2 



(10) 

Note that by stating 10 we do not keep the normaliza- 
tion of the Wigner function. Indeed, a finite fraction of 
the particles may hit the grating and may be removed 
from the flux. Therefore it is convenient to work with an 
unnormalized state and take care of the normalization 
only in the end. 

With the transformations JSJ) and we can proceed 
to describe the Talbot-Lau effect in a general framework. 



B. The Talbot-Lau setup 

In the Talbot-Lau setup a monochromatic beam passes 
three vertical gratings that are separated by the distances 
Li and Li. Since the particle stream is effectively uncolli- 
mated in front of the first grating its Wigner function for 
the transverse degrees of freedom is uniform. If we start 
with the (improper) normalization wo(r. p) = 1 then l|l(J|) 
yields the Wigner function after the first grating 



Wl (r,p) = |*i(r)| s 



(11) 



The free unitary evolution by a distance L\, followed by 
a passage through a grating (with convolution kernel T) 
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FIG. 1: The symmetric Talbot-Lau setup consists of three 
parallel gratings separated by equal distances L. Near-field 
interference of an uncollimated beam from the left may lead to 
a density pattern at the position of the third grating that can 
be observed by modulation with the lateral grating position 



and another free evolution by a distance L 2 , leads to the 
general expression 



w(r,p) = J dq 



t 1 [T-2-(L 1 +L 2 ) + —L l 

Pz Pz 



x T r 



Pz 



-L2 ,q 



(12) 



The particle density at position z — L\ + L 2 is obtained 
by integrating the momentum variable. It can be written 
as 

w{v)= y" W (r,p)dp = y"dr 1 |t 1 (r 1 )| 2 / l (r;r 1 ) (13) 
with 



""'■"[it 

Pz L x Li 



(14) 



As mentioned above the Talbot-Lau effect operates in 
the near-field regime, where the fact that the gratings 
have a finite lateral extension does not play a role; it 
only affects the overall count rate. It is therefore permis- 
sible to describe the gratings by idealized functions that 
are periodic on an infinite plane. Moreover, since the 
setup is invariant with respect to changes in the vertical 
position, it is sufficient to consider the Wigner function 
and the grating transmissions only as a function of the 
"horizontal" coordinates (r x ,p x ) = (x,p), see Fig. ^ 

We take the first grating to have a period d\, 
and its transmission function to be given by t\(x) = 
^2 dm exp(27rima;/di). Likewise, the second grating has 



the Fourier coefficients b m and the period d 2 
fore, 



d. There- 



sa;)! 2 = ^A^expf 2tt^ 



with A 



jez 



and 
T(x,p) 



(15) 



bjb*_ e exp 



hit 



2j-i 



. ' :Hi: 

In order to simplify the discussion and to avoid an 
overly complicated notation we focus on the symmetric 
Talbot-Lau setup, which is the most imp ortant one in 
practice; for the asymmetric setup see |2^| and note [70| . 
In this case the gratings are set at an equal distance L\ — 
L 2 = L and the periods di and d of the first and second 
grating are related by d\ = 2d/r, with r 6 N. (The 
case of equal grating periods, r = 2, is most common 
dill El El ) The state ^ now reads 



w(x,p) =- ^2 A i b j b j 



x exp^27ri 
x exp ( m£(2j 



d 
r L 
] 2T X 



Here we introduced the Talbot length 



x . . sLp 
m|— — 2ir\(r£ + m; 



d p z 



d 2 p z 
2-kK 



A 



(17) 



(18) 



in terms of the grating period d and the de Broglie wave- 
length A = h/p z . The Talbot length is the proper scale to 
distinguish near-field (Fresnel) diffraction from far-field 
(Fraunhofer) diffraction for a given wavelength. It gives 
the distance behind a grating where the diffraction peaks 
of a collimated, passing beam have a lateral separation 
equal to the grating period d. 

To get the particle density in the beam from (|17|l we 
integrate over the momentum p, which picks up the m = 
— rl components |7l| . 



w(x) = J dpw(x,p) 



oc 



eel. 



i A) exp( 27 ri£| 



with Fourier components [2c 



B W 



bib*_ m exp 



2jm L 



Ly 



(19) 



(20) 



Equation (|19f) predicts a density pattern which has the 
same period d\ as the first grating. Often the spatial 
resolution of detectors is too poor to detect these density 
oscillations directly in an experiment. However, an indi- 
rect observation is possible with the help of a third grat- 
ing with period d\. If put at the position of the density 
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pattern it modulates the total transmission as a function 
of its lateral position x s . The integrated transmission, 
which is much easier to detect, is then given by 



S{x s ) 



dp dq dx w(x,p — q)T(x 
dx w(x) | ^3 (x — x s )\ 2 . 



(21) 



If we choose the first and third grating to be identical, 
£3(2;) = t\(x), the expected periodic signal is given by 
the expression 



S{x s ) cx Y,( A *i) 2B er ) ^p( 2irie T 1 



(22) 



For symmetric gratings this modulation signal has a vis- 
ibility m 



Vc 



S n 



A 2n- 

n=l 



-\ D 2nr- 



Sn 



2 /i £i 



(23) 



n=l 



2nr 



which serves as the prime characterization of the inter- 
ference pattern. 

It is clear from the definition of the coefficients 
that the interference pattern Ijl9|l and the visibility l|23|) 
depend strongly on both the wavelength A and on the 
separation L between the gratings. As evident from l|2(J|) 
it is indeed the product of the two quantities which de- 
termines the pattern, since L/L\ = LA/d 2 . 

However, the detection of a periodic signal alone does 
not prove necessarily that quantum interference occurred 
in the experiment because a certain density pattern may 
also be expected from a generalized Moire effect. To es- 
tablish the observation of quantum interference one must 
show that the observed visibility differs significantly from 
the classical expectation. It is therefore important to 
have a reliable quantitative prediction for the classical 
expectation as well. 



C. The classical expectation 

With the results for the Wigner function at hand it is 
straightforward to repeat the calculation using classical 
phase space dynamics. The classical phase space density 
/(r,p) transforms under free evolution like the Wigner 
function JSJ according to 



/i(r,p) = /o ( r- — p,p 

Pz 



(24) 



In contrast to I|1U|) . the convolution kernel for passing 
through an ideal amplitude grating is now given by 



Tci(r,p) 



(2tt/i) 2 



dAe ipA / ft |f(r)| 2 = |t(r)| 2 <5(p) 



(25) 



which leads to 



/'(r,p) = |t(r)| 2 /(r, P ) 



(26) 



At a distance L 2 after the second grating this yields a 
phase space distribution 



/(r,p) 



h [r-—{L 1 + L 2 

Pz 



t 2 \r L 2 

Pz 



(27) 

which can also be obtained from the quantum result by 
replacing T by T c \ in Ijl2|l . It follows that the classical 
density pattern in front of the third grating is given by 



fix) 



with 



J dpf(x, P ) cx ^AlB^exp^iQ (28) 



B (o) 



5>. 



■ b* 



(29) 



The comparison with 119|) shows that the quantum and 
the classical results have the same form, but differ in the 
Fourier components B m . Of course, the classical Fourier 
components do not depend on the de Broglie wavelength. 
Nonetheless, the flS' may be viewed as the short-wave 
limit L/L\ — > of the quantum Fourier coefficients Bm\ 
which is already indicated by the notation. 

It follows immediately that the classical prediction for 
the signal is obtained from (|22|l by replacing the wave- 
length dependent Fourier components Bm by the Bm , 



S cl ( Xs ) CX J2( A i) 2B ir e M 27T[e T 1 



(30) 



Likewise one can show that the visibility of the classical 
signal is given, for symmetric gratings, by 



LXJ 

V A 2 R (0) 



AlB^ + £ A\ n B 



(o) 

2nr 



(31) 



Let us focus on the important case of equal grating pe- 
riods for a moment, i.e., r = 2. In this case only the even 
Fourier components B 2m are needed. If the separation L 
between the gratings is set to an integer multiple of the 
Talbot length, then it is easy to convince oneself that for 
any ideal grating t 2 (x) 



B 2 m - ° 



(0) 

2m 



if L/L x e N 



(32) 



that is, the quantum and the classical evolution yield 
identical predictions for the density pattern and the ob- 
served visibility. This shows clearly that the observation 
of the integer Talbot-Lau effect alone does not prove the 
wave nature of the beam particles. 
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Open fraction f 



FIG. 2: Talbot-Lau visibilities for ideal binary gratings as a 
function of the open fraction, (a) L — L\ (b) L — 0.9L\ 
(c) L = 0.8L\. The corresponding classical visibilities are all 
identical with curve (a). 



However, unlike their classical counterparts, the quan- 
tum Fourier components display a strong wavelength de- 
pendence. Therefore, distinctively different results are 
obtained in the classical and the quantum calculation 
for separations which differ from the integer Talbot cri- 
terion, L/Lx N, or equivalently for detuned particle 
wavelengths, A 7^ nd 2 /L,n € N. This can be seen in 
Fig. |21 where we show the quantum and classical visibil- 
ities for identical, ideal binary gratings as a function of 
their open fraction / (the ratio of slit width to grating 
period). 

As predicted by \.'>'2\) the quantum and classical results 
are identical for L/L\ = 1 and given by Fig. Ufa). For 
L/L x = 0.9 (Fig.GJb)) and L/L x = 0.8 (Fig.^c)) on the 
other hand, the Talbot-Lau visibilities differ markedly 
while the classical predictions remain on Fig. Ufa) . The 
distinction between the classical and the quantum pre- 
dictions is most pronounced for an opening fraction of 
0.5, where the classical contrast vanishes. The quantum 
calculation yields significant visibilities for these gratings, 
of 14.7% at L/L x = 0.9 and 25.4% at L/Lx = 0.8, re- 
spectively. 

D. Finite longitudinal coherence 

In the above calculations the particle beam that enters 
the interferometer was assumed to have a fixed velocity 
in the z-direction and to be completely uncollimated in 
the transverse direction. Of course this is an idealization 
that is in many respects as unrealistic as the familiar 
assumption of a perfectly coherent plane wave. Realistic 
particle beams are characterized by a finite longitudinal 



coherence and show some correlations in the transverse 
direction. 

The particle beams used in matter wave interferometry 
are usually generated by an effusive or supersonic expan- 
sion into a vacuum chamber |27| . By means of additional 
skimmers or collimators the beam is restricted to a well- 
defined "longitudinal" direction. The beam is stationary, 
and as a consequence the longitudinal momenta show no 
off-diagonal elements in their density matrix |28| . They 
are completely characterized by the longitudinal momen- 
tum distribution g(p z ) p9l |. 

The transverse momenta are much smaller than the 
longitudinal ones and can be taken to be uncorrelated 
with the longitudinal velocity. The transverse coherence 
is determined by the source aperture |30|. and it could 
be calculated with the van Cittert-Zernike theorem [3(j 
if the aperture size were small compared to the length 
scale in question. However, in the Talbot-Lau setup the 
source aperture is much larger than the spacing between 
the grating slits. As a consequence, diffraction at the 
first grating cannot be observed and it is permissible to 
approximate the transverse degrees of freedom as com- 
pletely uncollimated. In front of the first grating the bulk 
of the beam is therefore appropriately characterized by 
the Wigner function, 

VF bC am(R,P) = .g(P-e z ), (33) 

for R = r + ze z with z < 0. This description does not 
account for the edges of the beam and the cut-off at larger 
transverse momenta, which is why it cannot be properly 
normalized. Fortunately, it is not necessary to include 
the full beam profile in the treatment of the Talbot-Lau 
effect. As shown below only the interference of paths 
through a few neighboring slits is relevant for the effect, 
so that the transverse variation of the total current can 
be neglected. 

Formally, the beam ilMMI) can be written as a convex 
sum of states w(r, p) = 1 that are uniform in the trans- 
verse coordinates r and p, and that have a fixed longitu- 
dinal momentum p z . 

M / boam(r+ze 2 , p+p 2 e 2 ) = J 5(p z -p' z ) w(j^p^ g{p' z ) dp' z 

(34) 

Those are the states w(r, p) = 1 that we started out with 
in Sect. Ill 431 Since a sequence of grating transmissions 
© and free evolutions JHJ does not affect the dependence 
on p z the general stationary state at z' > is given by 

W / bcam(r + z'e z , p + p z e z ) = w(z'; r, p)g (p z ) (35) 

and the transverse position density reads 

w(r) = J dpdp z W hcam (r + z'e z ,p+p z e z ) 

= J dpdp z w{z'; r, p)g (p z ) . (36) 
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It follows that the finite longitudinal coherence in the 
beam is completely accounted for by averaging the results 
for a fixed velocity derived in Sects. Ill Bl and III CI over 
the longitudinal velocity distribution. In particular, the 
modulation signals l|22(l and (|30|l are given by 



(S(x s )) 



dp* g(pz)S(x s 



(37) 



If the detection signal is proportional to the flux the av- 
erage involves the longitudinal velocity component, 



(S(x B )) = 



JcT A Pz g(p z )p z S(x s 
J °° dp z p z g(p z ) 



(38) 



Note that in either case the visibilities are not obtained 
by a simple average of (12,'iH and (|31|l . but have to be 
calculated from the averaged signal. 



Here and below the tilde is used to indicate quantities 
which have an additional velocity dependence due to the 
grating interaction. Accordingly, for non-ideal gratings 
the convolution kernel l|l(JI) is replaced by 



T (x,p) = I dqT(x,p~ q)T v (x,q) 



with 

T v (x,q) 



x exp — i 



2irh 
. rripb 

Pzfr 



dAe iqA/h 



(40) 



(41) 



It follows that the quantum expressions (|19fl for the den- 
sity pattern and (|23|l for the signal visibility still hold 
after the replacement 



III. THE INFLUENCE OF REALISTIC 
GRATINGS 

So far it has been assumed that the gratings are ideal in 
the sense that their thickness could be neglected. How- 
ever, real gratings have a finite thickness and the time 
of interaction between the particle and the grating de- 
pends on the velocity v z — p z /m p of the beam parti- 
cles. This introduces a velocity dependence of the grat- 
ing function both in the classical and in the quantum 
treatments. Generally speaking, the Talbot-Lau effect 
is affected more strongly by the grating forces than far- 
field diffraction j^l), since the near- field interference is 
characterized by smaller phase shifts. This was seen 
in recent experiments with beams of large molecules 



B m ] - 5&> = £ y/ ,„ exp L?l^£!tjL\ (42) 

with the modified Fourier components b m — ^2 bjCm-j 
and 

c m = \ J\ e— /d ex P (-i^^) dx . (43) 

As mentioned above the presence of V(x) introduces a 
velocity dependence also on the classical level. This can 
be seen by considering the local approximation of l|41() 
where terms of the order d% V(x)A 3 are neglected in the 
exponent, 



A. The grating interaction 

In order to account for the effect of a finite grating 
thickness b we consider an additional interaction poten- 
tial V(x) that acts while the particle is traversing the 
grating. In order to avoid a more detailed z-dependence 
in the potential we average over the the surface roughness 
and assume that the grating walls are parallel to the opti- 
cal axis and that edge effects can be neglected. In the case 
of tilted walls one can introduce an effective slit width, as 
discussed in . Moreover, it is known [33103 that grat- 
ing interaction effects are usually well described by the 
eikonal approximation. There, the additional quantum 
phase due to the interaction potential V(x) is obtained 
by integrating the action along a straight path. Accord- 
ingly, if we take the binary function t{x) to describe the 
material grating the complete grating function is given 
by 



, .m D bV(x) 
t{x) = i(a;)exp -i^ ^ 



(39) 



rripb d 
Pz dx 



V(x) 



h. 



= S [q + — — — V(x y 
Pz dx 



(44) 



According to l|4U|) this yields a classical convolution ker- 
nel 



f cl (x,p) = \t 2 (x)\ 2 6(p+ 1 ^^-V(x) 
V Pz da; 

= \t 2 (x)\ 2 5(p-Q(x)) 



(45) 



which indicates that the eikonal approximation corre- 
sponds on the classical level to the momentum kick 
Q(x) = — d x V(x) x b/v z obtained by multiplying the 
constant classical force at a fixed position x with the 
interaction time. Accordingly, the classical phase space 
distribution changes as f'{x,p) — \t(x)\ 2 f(x,p — Q(x)) 
when passing a grating. Using this transformation and 
the periodicity of Q one finds that the classical expres- 
sions for the density pattern and for the signal visibility 
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assume the forms (|28ll and (|31|l as in the ideal case. One 
merely has to replace the Fourier components by 



with 



s-iTl 



e -27rimx/d exp 



LQ(x) 



(46) 



dx . (47) 



We note that the modifications of the Fourier components 
given by and (f4*6*j) describe the quantum and the clas- 
sical interactions on the same degree of approximation. 
Clearly, the fact that the interaction with the grating 
is treated equally in the quantum and in the classical 
descriptions is an important requirement for identifying 
quantum interference in an experimental observation. 

Before turning to realistic descriptions for the grating 
interaction we note that it is in general not necessary 
to include the grating interaction at the first and third 
gratings. This is clear from the fact that the Talbot- 
Lau setup is sensitive to diffraction only at the second 
grating, while the others merely serve to modulate the 
flux. Formally, it can be seen from the expression for the 
observed signal, Equation (|21[1 with 1|13|) . It depends only 
on the squared moduli |ii(a;)| 2 and ^(x)! 2 of the first 
and third grating function, which are not affected by the 
phase shift in (|39[1 . Only for very strong potentials, where 
(|39[1 is no longer valid, may the interaction effectively 
reduce the slit width and thus become relevant to the 
first and third grating. 



B. Material gratings 

A neutral particle will in general experience an attrac- 
tive van der Waals force if placed in the vicinity of a 
surface. A simple, but quite realistic description for non- 
polar quantum objects is given by the static London dis- 
persion force which acts between a quantum object and 
a flat wall. If A is the distance to the wall it give s rise to 
the potential U(A) = -C 3 /A 3 , with C 3 > 0^|. 

For simple wall materials the interaction constant C3 
can be found in the literature for many atoms and a 
number of small molecule s 1361. I n ge neral it is obtained 
from the Lifshitz formula 35, 37, 381 



h r°° ^M-i , 

C 3 = — / a (iuj)— — — du 
4tt J e(iw) + 1 



(48) 



by using either experimental data (e.g. absorption spec- 
tra) or appropriate models for the dynamic polarizability 
a of the particle j?^ and for the bulk dielectric function 
e of the grating material, respectively. Often Prude- t ype 
models for a and e are considered sufficient. In [33l l34j] 
the interaction with material gratings was studied in an 
interference experiment and found to be in good agree- 
ment with the assumption of a London dispersion force. 



However, at large distances retardation effects may be- 
come important. They show up if the separation A be- 
tween the particle and the grating wall is comparable 
to the wavelength corresponding to those virtual transi- 
tions in the particle that contribute with a large oscillator 
strength. In the case of an ideal metal the potential is 
described by the Casimir-Polder formula 39]. For large 
distances it has the asymptotic form U(A) — — C4/A 4 
with the constant 



Ca 



3hc 
87 



0(0) 



(49) 



given by the static polarizability a(0) of the particle. The 
case of more realistic grating materials and arbitrary dis- 
tances is covered by the theory of Wylie and Sipe |40ll4l| . 
It shows that in the case of real metals the asymptotic 
form of the interaction potential does not depend on the 
metal and is identical to the ideal case. For dielectrics 
the limiting form depends also on the fourth power of the 
distance, but with a reduced interaction constant, 



Ca — Ca 



(1 + 2u 2 )r p (u) - r s (u) u 



(l + u 2)5/2 



-du. (50) 



The reduction depends on the static dielectric constant 
e(0) via the Fresnel coefficients 



r P (u) 



and 



r s (u) 



vMO) 



(51) 



(52) 



Figure |21 gives the value of the dielectric reduction factor 
for 1 < e(0) < 100. 

Whether the exact position dependence of the retarded 
force must be used depends on the physical situation in 
the particular interferometric setup. In most experiments 
realized so far it was sufficient to use either the static van 
der Waals interaction Ij48(l 31] or the long range limit of 
the Casimir-Polder force, Equation l|49l) . 

Figure 0] shows the typical effect of a finite grating 
interaction, and should be compared to the results for 
the ideal grating in Fig. [21 Here we assume a particle 
with mass 1000 amu, a van der Waals interaction with 
C3 = lOmeVnm 6 , and we take gratings with a period 
of d = 1 /xm and a thickness of b — 0.2 /im separated 
by a distance of L = 0.2 m. One observes that the ex- 
pected quantum visibilities deviate noticeably from the 
ideal case. Moreover, the classical expectations (given 
by the dashed lines) differ completely from the ideal ex- 
pectation and display now a weak velocity dependence. 
At an open fraction of 0.5 they now yield a finite con- 
trast amounting to 12.4% — 12.9% for the three settings. 
The respective quantum expectations are also larger than 
the in the force-free calculation. (They increased from 
14.7% to 32.9% at L/L x = 0.9 and from 25.4% to 35% 
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FIG. 3: Reduction factor C4/C4, see Eq. H5UII . of the long- 
range interaction for dielectric gratings as compared to ideal 
metals (semi- logarithmic scale). For typical grating materials 
the static dielectric constant e (0) is less than 4, leading to a 
reduction below one half. 
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FIG. 4: Talbot-Lau visibilities for gratings with a van der 
Waals interaction as a function of the open fraction, (a) L = 
L\ (b) L = 0.9L\ (c) L — 0.8L\. The corresponding classical 
visibilities are given by the curves (d), (e), (f), respectively. 



at L/L\ = 0.8.) This is a typical phenomenon. The at- 
tractive force tends to act as if the open fraction of the 
grating was decreased. 



C. Gratings of light 

It is clear from Equations l|13(l and (|21|l that the first 
and the third gratings in the Talbot-Lau setup must be 
absorptive to generate an observable contrast pattern. 
However, as discussed in [23j the second grating may be 
a pure phase grating as well. Such a mixed interferome- 
ter can be realized by the off-resonant interaction with a 
standing light wave p2|> see E3 f° r the laser diffraction 
of large molecules. 

If we take a TEMoo mode of wavelength Al = 27r/fcL 
and waist w produced by a laser of power Pl then the 
dipole force leads to the phase shift 



t{x) = exp i 



j7r ^ LQ; " cos 2 (k^x) 
hcv z w 



(53) 



as follows from an integration over the gaussian beam at 
central passage. Here a u is the scalar polarizability [73] 
of the particle at Laser frequency cjl = c&l, where we 
assume that Im(a UJ ) — so that photon absorption can 
be neglected. The classical momentum kick in (|45H that 
corresponds to the dipole force reads 

Q(x) = v^ 8 ^"" kh sin{2k L x) . (54) 

CVyW 

Using these expressions in (I39II and (|47ll one obtains the 
predictions for the quantum and classical density pat- 
terns produced by a standing light wave in the same way 
as with material gratings. 

In the present section we showed that the Wigner de- 
scription of the Talbot-Lau effect permits the effects of 
the grating interaction to be incorporated easily, in terms 
of a simple modification of the Fourier coefficients. As 
discussed in the next section the effect of decoherence 
can be similarly incorporated into the formalism. 



IV. ACCOUNTING FOR DECOHERENCE 

Having formulated the Talbot-Lau dynamics in the 
Wigner representation it is now easy to include the effects 
of decoherence. More specifically, we consider the marko- 
vian interaction of the interfering particle with other, un- 
observed degrees of freedom (the environment) fiol Hl|. 
The resulting formation of quantum correlations (or en- 
tanglement) between the particle and the environment 
leads to a loss of coherence in the particle state that may 
be understood from the fact that a measurement of the 
environmental degree of freedom could reveal (partial) 
which-way information on the particle's whereabouts. 

We note that a number of studies have been under- 
taken recently that describe a loss of visibility in mat- 
ter wave interference [Iii liol |4^ | . Here we focus on the 
Talbot-Lau effect and on a formulation that is sufficiently 
realistic to permit quantitative predictions about exper- 
iments with mesoscopic bodies 



preciicti 

(am. 
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Two important decoherence mechanisms for large, in- 
terfering particles are collisions with background gas par- 
ticles and the thermal emission of electromagnetic radia- 
tion. Both effects may be treated in the Markov approx- 
imation, which implies that the effect of the environmen- 
tal coupling can be described by independent, separate 
events (such as the emission of a photon or the collision 
with a gas particle). 



magnitude of p z . Our basic approximation is now to as- 
sume that the width of 77 is small compared to the scale 
over which Wbeam varies in p z . This assumption is par- 
ticularly unproblcmatic in the Talbot-Lau setup, where 
the sensitivity to changes in the longitudinal momentum 
is rather weak. It follows that the new state of the trans- 
verse coordinates is approximately given by integrating 
the full state with respect to the longitudinal momentum, 



A. The effect of a single decoherence event w'(z', r, p) ~ J dp' z W^ cam (r + ze z , p + p' z e z ) . (60) 



The change in the state of the interfering particle due 
to a single event can be obtained by performing a partial 
trace over the entangled state with respect to the un- 
observed degrees of freedom. For particles with a large 
mass and for the decoherence mechanisms considered in 
this article the density matrix in position representation 
changes just by a multiplication, 



e'(Ri,R 2 ) = e(Ri,R 2 )?7(Ri 



R 2 



(55) 



The factor ?j(Ri — R2), which may be called the deco- 
herence function, describes the decay of the off-diagonal 
elements (the coherences) of g due to a single event. In 
Sect. El we will derive the form of realistic decoherence 
functions for the most important decoherence mecha- 
nisms. For the time being it is sufficient to note that 
the conservation of the trace in 15511 ensures that 



lim n(Ri 

Ri^R 2 



R 2 



1, 



(56) 



so that the diagonal elements of the state are unchanged 
by H55JI. Moreover, the hermiticity of g implies 77 (— R) = 
77* (R) , and from the fact that the purity cannot be in- 
creased by a partial trace it follows that |r/(Ri — R 2 )| < 1. 

If the state is expressed in terms of the Wigner function 
its change reads 



W'(K,P) = J dQf?(Q)W(R,P-Q). 



(57) 



with 7j(Q) the Fourier transform of the decoherence 
function, 



-Ur / dR e- iQR / fi 7j (R) 
irh) 3 J 



(2nh) 



(58) 



Clearly, the effect of a decoherence event on the Wigner 
function is to smear it out in the momentum direction. 

As discussed in Sect. Ill Dl the coherently evolving, sta- 
tionary state of the beam in a Talbot-Lau interferometer 
is described by the function 

W hcam (r + ze z , p + p z e z ) = w(z; r, p)g(p z ). (59) 

In a typical setup the grating constant and the grating 
separation differ by six orders of magnitude so that p 
varies on a scale in <|59[) that is much smaller than the 



Inserting l|57() and (|59|l yields the change in the state of 
the transverse coordinates 



w'(z; r, p) = J dq 77 2d (q) w (z; r, p - q) (61) 



with 



»72d(q) = J d<? 2 77 (q 4- q z e z 



(62) 



It follows from (1621 that in the position representation of 
the transverse state, /j(ri,r 2 ), the decoherence function 
enters without modification. 



p'(r 1 ,T 2 ) = p(ri,r 2 ) J dq 



,iq(ri-r 2 )/fi 



I2d(q) 



p(ri,r 2 ) 77(1-1 -r 2 ), 



(63) 



but restricted to the xy-plane, 77 (r) = 77 (r + 0e z ) . 

Using the Wigner function it is now possible to eval- 
uate the effect of a single decoherence event that takes 
place at a distance z behind the first grating. One merely 
propagates the state to the longitudinal position z using 
the coherent transformations JHJ and then applies 
(|61|l . and propagates the state over the remaining dis- 
tance to the third grating. Within the paraxial approx- 
imation no additional error is introduced by this pro- 
cedure since the composition property of the evolution 
holds exactly. The result takes a simple form (both for 
z < L and for z > L) once the momentum is integrated 
to yield the position density w (r) in front of the third 
grating. It is given by an integral of the form l|13fl . 



w (r) 



dri|fi (n)| a ft,(r;n; 



(64) 



where the coherent kernel /i(r,ri) from l|14|) is replaced 
by 



/i z (r;ri) = / dq77 2d (q) h[r,ri - 



L\ 



Pz 

(65) 

This shows clearly that close to the first and to the third 
grating, at z — and z — 2L, a decoherence event will 
not affect the interference pattern while, for monotoni- 
cally decreasing 77 2( j(q), the interference is most strongly 
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affected by decoherence events that take place in the 
vicinity of the second grating, at z = L. This is consistent 
with the notion that in the Talbot-Lau setup diffraction 
takes place only at the second grating, while the first 
grating acts as an array of coherence slits. 

We take the grating function again to be periodic in x 
(with period d) and uniform in y. This means that the dis- 
cussion can be confined to the rr-coordinate like in Sect. 
Ill Bl Using the Fourier decomposition 1|16|) one finds that 
the coherent kernel (|14fl reads in the one-dimensional case 

h(x;xx) = J dydy 1 h(xe x + ye y ]x 1 e x + y 1 e y ) (66) 



Pz 

2Lh 



exp | 27ri//f 



x + Xi\ ~ m 
2d 



It follows with H65[l that in the presence of a decoherence 
event the kernel takes the form 



I i \ ft V ■ x + x i 

h z {x;x x ) =— 2^ CX P [ 

d L — \z — L 



2Lh 

771 

x B^r, I -/,/- 



27rim- 
V 2d 



ex , (67) 



where the three dimensional decoherence function from 
(|55|l enters with its dependence along the x-axis. The 
comparison with (|66[) shows that the modified interfer- 
ence pattern corresponding to a single decoherence event 



at position z is completely described by a modification 
of the coherent Fourier components l|42l) 



dL- \z-L\ 



-to—- 
2 



(68) 



B. An alternative to the master equation 

One can now account for probabalistically occurring 
decoherence events by considering the change in the final 
interference pattern due to events that occur with rate 
R(z) in the interval (z; z + dz). It follows from (|68|l that 
the corresponding Fourier coefficients satisfy the differ- 
ential equation 



±BW 

dz m 



R(z) 



dL- \z-L\ 



-TO— - 

2 



_ (69) 

It describes the change of the interference pattern with an 
increasing size of the interval where decoherence events 
may occur. The integration of (|69|l over the whole range 
z € (0; 2L) of admitted decoherence then yields the co- 
efficients characterizing the modified pattern. They are 
given by 







-f R,3) 


1 — */ ( 







dL-\z-L\ 



dz 



(70) 



r 



with fim' the coefficients of the coherent evolution. This 
is the central result of this section. It shows that the 
effects of markovian decoherence of the form (I55|l can be 
calculated analytically if the setup is insensitive to longi- 
tudinal correlations as in the Talbot-Lau interferometer. 
It follows immediately that the position density and the 
visibility of the modulation signal are given by the formu- 
las H19(l and Ij23(l , respectively, if the coherent coefficients 
Bm' are replaced by those of the incoherent evolution 
GDI- 

The result l|70|l can be easily generalized to the asym- 
metric Talbot-Lau interferometer. The case of several 
independent decoherence mechanisms is also easily in- 
corporated. The resulting interference pattern is then 
characterized by a product of the corresponding expo- 
nentials in (|70p. 

It is important to note that the basic Fourier com- 
ponents to = are not affected by decoherence, since 
77 (0) = 1. This shows that the mean count rate does not 
change due to the presence of decoherence, as is to be 



expected from the conservation of the norm in (|55|l . The 
reduction of the observed visibility assumes a compact 
form if the modulation signal (1221 is (approximately) si- 
nusoidal, as is typically the case for gratings with an open 
fraction of / ~ 0.5. Then only the coefficients -Bg^and 

B^ contribute to the visibility if the grating periods are 
equal, r = 2. With 77 (0) = 1 it follows that the reduced 
visibility is given by 



V = V exp - 



2L 



R{z) 



1-7? -d 



L - 



L\ 



(71) 

where Vo indicates the visibility in the absence of de- 
coherence. This formula is particularly intuitive if the 
Talbot criterion is met, L = L\. Then the argument of rj 
contains the separation of two paths that start and end 
at common points and pass the second grating through 
neighboring slits. At z — L it is equal to the grating con- 
stant d, which shows that the Talbot-Lau interference 
with equal gratings is based on the interference through 
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neighboring slits. Also at other positions z a reduced 
magnitude of 77 suppresses the visibility whenever the 
change in the environmental state is able to resolve the 
corresponding path separation. Higher orders of the Tal- 
bot effect, L = mL\ with m € N, correspond to multiple 
slit separations md. For longitudinal velocities that de- 
viate from the Talbot criterion with L ^ rnL\ the argu- 
ment is replaced by an "effective" path separation. 

It should be emphasized that our derivation of i|70|) 
and l|71f) is rather different from solving the markovian 
master equation corresponding to the decoherence mech- 
anism. In Sect. IVII we obtain a solution of the master 
equation corresponding to decoherence of the type l|55|) 
for a general interfering state in the paraxial approxi- 
mation. An expression analogous to (|71(l is found there, 
albeit in a time-dependent formulation, see i|12U|) . This 
vindicates our approximation (|60|l . 

The present formulation has the particular advantage 
that the rate R of decoherence events and their effect 
77 appear separately in the equation. This might seem 
to be a complication, since these two quantities must be 
calculated independently by quantum mechanical means. 
However, they are often needed with different degrees of 
sophistication. For example, often one must take into 
account the position dependence of the rate. This is eas- 
ily incorporated in the present framework, while solving 
a corresponding master equation would be incomparably 
more complicated. 

C. Quantum decoherence vs. a classical stochastic 
process 

Having treated the effect of environmental coupling on 
the quantum evolution, we can now turn to its effects 
in the classical description. In Sect. Ill CI the classical 
expectation was calculated in terms of the phase space 
density /(r,p). The close analogy between the quantum 
problem and the classical calculation allows to map the 
Wigner representation of a decoherence event Q61[) to the 
classical description. It follows from (|61[) that the effect 
of a decoherence event can be interpreted on the classical 
level as a probabilistic momentum kick, 

/'(r, p) = J dq fj 2d (q) / (r, p - q) . (72) 

Indeed, the properties 77 (0) = 1 and 77 (— r) = 77* (r) of 
the decoherence function imply that 7/2d has the features 



of a probability density, ?72d(q) > and / dq ?y"2d (q) = 1- 
From the close analogy of the classical and the quantum 
expressions for the free evolution and the passage through 
a grating it is easy to see that the above derivation of the 
modified pattern holds in the classical formulation as well 
if one replaces the quantum coefficients B$ in l|7U|l by 
their classical counterparts Bm ■ 

From this one might be led to conclude that the deco- 
herence described in (|55|) was a "classical effect" . In our 
view this would be a misinterpretation, since a probabilis- 
tic formulation is possible only if the Wigner function is 
non-negative everywhere, that is, if it cannot be distin- 
guished from a classical probability distribution. If the 
Wigner function is negative in some parts, as is the case 
for an interfering state, any stochastic interpretation is 
invalidated by the occurring flux of a "negative probabil- 
ity" . Notwithstanding this, once the motional state has 
turned into a classical state without negativities in the 
Wigner function the additional loss of visibility in the 
quantum description is indeed indistinguishable from a 
corresponding classical stochastic process. 



V. REALISTIC DECOHERENCE FUNCTIONS 

In the following we discuss the form of realistic deco- 
herence functions that can be used to obtain quantitative 
predictions on the effects of decoherence in matter wave 
experiments. We focus on the most important mecha- 
nisms for large, massive objects, namely, collisions with 
particles from the background gas and the emission of 
heat radiation. We note that simple estimates of these 
effects on material particles can be found in 0, HE E3 • 



A. Decoherence by collisions with gas particles 

A very important source of decoherence is the unavoid- 
able presence of a background gas in the experimental 
apparatus. Typically, the mass of the interfering parti- 
cles is much larger than the mass of the gas particles, 
m p 3> m g , and the interaction is of the monopole type. 
In this case the decoherence function reads 1471 150, 51 1 



77(Ri,R 2 ) = tr gas jexp f-iP gas ^J §1 cxp ^iP gas So cxp ^-iP gas ^-^ £ gas j (73) 

where P gas is the momentum operator of the gas particles and So the center-of-mass scattering operator. The trace 
over the scattered gas particle in l|73|l can be evaluated if it is in a (thermal) state that is diagonal in momentum and 
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characterized by the distribution /i gas (P) . One obtains [50ll5l| 



tfRi.Ba) = JdP Mgas (P) [l-/dP' (i_ e i(P-P')(^-^)M) IMl!| (P '|f |P>| 2 



(74) 



with To = i(l — So). Awkwardly, the last two expressions 
in (|74fl . which are indicated by JTq, involve two quan- 
tities that are arbitrarily large. One is the "quantiza- 
tion volume" f2, which originates from the normalization 
of the thermal state g gas , and the other is the square 
of the delta function appearing in the matrix element 
<P'|f |P) = f(P', P)S (P' - P) /(2nhP). Here /(P', P) 
is the scattering amplitude (which must not be confused 
with the classical phase space density from Sect. Ill Cll . 
Since the decoherence function is well defined by Equa- 
tion (|73|) these two infinite quantities must cancel if the 
limit H — > oo is taken properly. As argued in |5lj| physical 
consistency requirements lead to 



lim Xn = 



|/(P',P)| 2 5(P' -P) 



a(P) 



P 2 



(75) 



(76) 



(77) 



where a (P) is the total scattering cross section 

a(P) = Jdn\f(Pn,P)\ 2 . 
With the replacement l|75[) one gets 
^(Rx,R 2 ) = JdP^ 1 

X Jdn |/(COS(0))| 2 cKP-PnXRx-R*)^ 

As already anticipated in Equation (|55|l this function de- 
pends only on the position difference Ri — R2. For an 
isotropic distribution of the gas momenta, /i gas (P) — 
^gas (P) I (47rP 2 ) , the expression can be further simpli- 
fied noting that it depends only on the distance AR = 
IRi — R2I. One obtains 



V (AR) 



°° ir ^s(P) 

°(P) 



(78) 



x JdQ |/ (cos(6>))| 2 sinc^sin^ 



0\ 2PAR\ 



with sine (x) = sin (x) /x. The argument of the sine func- 
tion is equal to the momentum transfer during the colli- 
sion times the distance AR in units of ft. This indicates 
that whenever the change in the state of the gas parti- 
cle suffices to resolve the distance AR the corresponding- 
coherences in the motional state will be suppressed. 

Let us turn to the second ingredient to the decoherence 
formula 171|) . the scattering rate R(z). It is usually ex- 
pressed in terms of an effective cross section, R (z) = 



n(z)a e g, with n (z) the number density of the back- 
ground gas. For a constant density, and again m p 3> m g , 
the effective cross section depends only on the velocity of 
the interfering particle. It is given by 



|P/m g 



Ceff («p) = / dP ^ gas (P) a (|P - TO g V p e z |) 

J vp 

(79) 

as follows from the derivation of the Boltzmann equation. 

The most prominent interaction encountered in molec- 
ular scattering is the van der Waals force between polar- 
izable molecules. At the typical velocities in matter wave 
interferometry the scattering depends only on the long- 
range part of the interaction potential, U (r) = —Ce/r 6 , 
which is characterized by a single interaction constant 
C% . The total cross section is then independent of mass 
and given by jS^ 



a (jn sV ) = 



r (2/5) sin (tt/5) 



~8~ ft^ 



2/5 



(80) 



The integration in (|79|) can be done assuming a thermal 
distribution of the gas particles. The exact expression is 
given by a confluent hypergeometric function, as shown 
recently by Vacchini |5_3j. Here we note the asymptotic 
form of the effective cross section Q79JI for small velocities 
of the interfering particle. It reads 



4nT 



(9/10) / 



3irC 6 
5 sin (tt/5) V 2ft 

N < I + - I — ) + O 



2/5 -3/5 




(81) 



1/2 

with Vg — (2k^T/m s ) the most probable velocity in 
the gas. 

In principle, the interaction constant is given by the 
Casimir-Polder expression 



C 6 = 



3ft 



/6 = — / duo a„ (iuj) a„ (iuj) 

TT 



(82) 



involving the frequency dependent polarizabilities of the 
two particles [73]. However, often only the static polariz- 
abilities are available for larger molecules. In this case a 
fairly accurate estimate can be obtained from the Slater- 
Kirkwood expression |54j 



Co 



a s (0)a p (0) 



2 V47T£ m e y/ag (0) /N e + y/a p (0) /N p 



(83) 
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where N g and N p are the number of valence electrons 
of the gas molecules and the interfering particle, respec- 
tively. 

Let us stress again that in the present treatment the 
effect of a single collision (|78|l and the rate (|81ll are cal- 
culated separately, which is particularly useful if the two 
are needed at different degrees of accuracy. This was the 
case in the recent experiments on collisional decoherence 
[Hi. l55j where the localization took place on a scale that 
is by orders of magnitude smaller than the path sepa- 
ration. Consequently, 77 could be replaced by a simple 
Kronnecker-like function in i|7U|) . while the finite veloc- 
ity of the interfering particle within the thermal gas had 
to be taken into account properly. Corresponding mas- 
ter equations, based on the microscopic description of 
the scattering proc ess, are a subject of current research 
El El El El El El Although some of those are 
sufficiently detailed to describe the emergence of an effec- 
tive scattering rate (I81|l . their application to a description 
of the experiment would have been considerably more 
complicated than in the present treatment. 



B. Decoherence by thermal emission of radiation 

A second decoherence mechanism that is common to 
all macroscopic objects is the emission of heat radiation. 
It starts to play a role in matter wave interference if one 
considers macro-molecules or mesoscopic particles. Due 
to the large number of internal degrees of freedom, a 
thermodynamic description of the distribution of the in- 
ternal energy is unavoidable. Moreover, their coupling 
to the electromagnetic field is quasi-continuous. In gen- 
eral, the thermally emitted photons will reveal (partial) 
which-way information on the whereabouts of the inter- 
fering particle and thus lead to decoherence. 

We assume that the emission is isotropic, and that the 
walls of the apparatus, which absorb an emitted photon, 
are located in the far field where the photon's spatial 
detection probability is given by its momentum distri- 
bution. The conservation of the total momentum then 



suffices to determine the transformation of the particle's 
density operator that would be obtained from a partial 
trace over the entangled state between photon and par- 
ticle. It follows that the change of the particle center of 
mass coordinate due to a single emissions given by 



dk !§^ (84) 



where Pk(k) is the probability distribution for photons 
with wave number k = |k| and C/k = exp(iRk) are the 
momentum translation operators. Note that it is not 
necessary to consider the change of the internal degrees 
of freedom of the particle, since their state does not get 
entangled with the center of mass; this would result only 
if the emission probability were position dependent. 

In position representation, £>(Ri,R.2) = (Ri|g|R2), 
the transformation l|84|) reduces the off-diagonal elements 
of the center-of-mass state 

g'(Ri,R a ) = e(R 1 ,R 2 )»7(Ri-R 2 ). (85) 

The corresponding decoherence function reads 



rj(AR) = —J^dX R\(X) sine (2,— 



(86) 



where the probability distribution pk was expressed in 
terms of the spectral photon emission rate, 



_ ... 2irR tot /2tt 
#a(A) = A2 Pk I y ] 



(87) 



and the total photon emission rate 

I* oc 

Rtot = / fl A (A)dA. 
Jo 



In the expression for the Fourier components (|70|l the 
total rate of decoherence events i?tot cancels because it 
gets multiplied by r\. One obtains 



2L 



dz 



dX Rx(X) 



SinC TO7T 



dL-\z-L\ 
A 



Ly 



This shows clearly how the fringe pattern gets blurred 
by heat radiation if it contains photons that have a suf- 
ficiently small wavelength to resolve the path separa- 
tion. The properties of the interfering particle enter only 
through the spectral emission rate R\ = R^ x |dw/dA|. 

For mesoscopic particles the spectral emission rate de- 
viates from the Planck law of a macroscopic black body 
for a number of reasons. First, the photon wavelengths 



are typically much larger than the radiating particle, 
which turns it into a colored emitter. The density of 
available transition matrix elements can be related to the 
absorption cross section [o(l | . Second, at internal energies 
where thermal emission is relevant the particle is usually 
not in thermal equilibrium with the radiation field, so 
that there is no induced emission. Third, the particle is 
not in contact with a heat bath, but the emission takes 
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place at a fixed internal energy E. Similarly to Einstein's 
derivation of the Planck law, these points lead to the ex- 
pression jH3 



Ru (w) duj 



7 cr a bs (E - Huj; uj) 



d(E-hu) 
d(E) 



duj. (89) 



The first term is proportional to the mode density. The 
mean oscillator strength is described by the absorp- 
tion cross section at frequency u> and internal energy 
E — fkj, and the ratio of the densities of state d(E) 
yields the statistical factor under a strong mixing as- 
sumption. The mean densities of states can be re- 
lated to the thermodynamic properties of the particle 
by a stationary phase evaluation of the inverse Laplace 
transform of its partition function [(J^, |6j| . This yields 
d (E) ~ exp (S (E) /fee) and therefore 



d(E-tiw) 
d{E) 



exp 



fku 
knT* 



1 



2C V 



(90) 



Here, the internal energy is conveniently expressed in 
terms of the micro-canonical temperature 



T* (E) = 



OS (E) 



dE 



(91) 



with S (E) the entropy. The value of T* is equal, up to 
small corrections, to that canonical temperature where 
the mean energy equals the internal energy. The second 
term in (|90[1 contains the heat capacity Cy of the particle. 
It is the leading correction due to the finite size of the 
internal heat bath. This term decreases with increasing 
size of the particle and l|90(l assumes the canonical form 
in the limit Cy — > oo. 

With equations (|89|l and H9U|) one is able to calcu- 
late the temperature dependent spectral emission rate 
R\(\,T) and the corresponding decoherence effect. At 
very small heat capacities the effect of cooling may have 
to be taken into account. It can be easily incorporated 
in the present framework through a position dependent 
temperature T (z) determined by the cooling formula 



4~T (z) = \— [ hujR u {uj,T)duj. 

dz v z C v J 



(92) 



We note that also scattering of photons may lead to 
decoherence, although room temperature photons will 
not limit matter wave interference in the foreseeable fu- 
ture. The deliberate scattering of a laser beam at inter- 
fering atoms at a resonant cross section was studied in 

[UlEllHEllfti 

Finally, we emphasize that all the calculations in this 
paper have been within the framework of conventional 
quantum mechanics. However, proposed extensions of 
that theory, which produce spontaneous localization of 
massive particles due to postulated "collapse" terms 
added to the Schrodinger equation, lead to an evolution 



of the density operator which mimics decoherence effects 
09 . Hence both the establishment of the framework we 
develop here, and accurate models for decoherence mech- 
anisms, are essential to ascertain whether or not any par- 
ticular proposed modification of conventional quantum 
mechanics can be ruled out by experimental data. We 
defer such applications to future articles with a focus on 
laboratory results. 



VI. EQUIVALENCE WITH THE MASTER 
EQUATION 

In this final section we show that the procedure to in- 
corporate decoherence that was used in Sect. lIVl is equiv- 
alent to solving the corresponding master equation in 
paraxial approximation. This is done by identifying the 
systematic corrections to the paraxial approximation in 
terms of the ratio between transverse and longitudinal 
momenta. 

Our starting point is the master equation for a free 
particle 



d 



1 



dt 8 ift 



2m, 



-,Q 



J dRdR' 7 (R — R') <?(R,R') |R)(R'| 



(93) 



with localization rate 7. It is valid in situations where 
the mass of the particle m p is sufficiently large so that 
the effect of the environmental coupling does not (yet) 
lead to thermalization. This equation is usually applica- 
ble in interferometric situations where one is interested 
in time scales that are much shorter than those of dissi- 
pation. In particular, it describes the effects of scattering 
of particles with a much smaller mass or the emission of 
photons. 

It follows from (|93|l that the corresponding Wigner 
function satisfies 

|- W (R, P; t) = - — V R H/ (R, P; t) 

at m p 

- [dP'*/(P') W(R,P-P';t) (94) 

with 7 the Fourier transform of the localization rate, 

7 (P) = -^-3 / dR e- iRP /" 7 (R) . (95) 
(27m) J 

Unlike in the previous sections, we describe the motion 
of the particle by a Wigner function that is properly nor- 
malized, 



J dRdP!F(R,P;<) = 1. 



(96) 
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A. Decoherence of an interfering state 

Now consider the usual scattering situation where 
the particle enters and leaves the grating region of 
an interferometer in a finite period of time, so that 
W(R, P;t — > ±00) = for all positions R of interest. 
It follows that 

/oo 
dt d t W (R, P;t) = W (R, P; +00) - W (R, P; -00) 

= 0. (97) 

As above, we take the z-axis as the longitudinal direction 
of the interferometer, 



R = r + ze z 
P = P+Pze z , 



(98) 



and denote the transverse positions and momenta by r = 
(x,y) and p = {p x ,p y ), respectively. 

At t = the particle is localized in the region z < 
and heading for the region z > where decoherence may 
occur, say, because there is a gas present. Moreover, we 
assume that at t = the particle is already in a non- 
classical motional state, for example because it has just 
passed a grating. The expected interference pattern is 
given by the position dependent detection probability in 
the z-plane which is obtained by integrating the longitu- 
dinal current density over time, 

Q(z;r) = dt dpdp z —W (r+ze z ,p+p z e z ;t) . 
J J m p 

("> 

In order to compare with the results from Sect. lIVl we are 
ultimately interested in the effect of decoherence on the 
Fourier transform of the interference pattern with respect 
to the transverse coordinates, 

Q(z;q)=-l- 5 /clre- i ^Q(z;r) 
(27m) J 



dp dp. 



Pz 



dre" iqr/fi 



x J dtW (r+ze z ,p+p z e z ;t) 
= J dpdp z S a (z,p z ;q,p) . 



(100) 



Here we introduced the auxiliary function S a . In 
order to obtain a differential equation for S a apply 
Jdr exp (-iqr/fi) Jdt[-] to Using and inte- 

grating by parts one finds 

— S a (z,p z ;q,p) = -i — ^-S a {z,p z ;q,p) (101) 

-m p f dp' z dp' ^ P + Pz * z S a (z, p z -p' z ;q,p- p') ■ 

J Pz- p' z 



In the absence of decoherence this differential equation 
is immediately integrated, 

S a (z,p x ;q,p) =exp (- iz< ^jf) s a (0,p*;q,p) (102) 
(for 7 = 0). 

This decoherence-free solution is used below to obtain a 
systematic approximation in the presence of decoherence. 
But first we introduce the Fourier transform of S a with 
respect to the longitudinal momentum referenced by a 
fixed characteristic momentum p z , 

Sb (z, C;q, p) := J dps exp (Kps/K) S a (z,p z +Pi;q,p) 

(103) 

so that 



Q (z;q) = dpS b (z,0;q,p) 



(104) 



The motivation for this definition is that we will assume 
that S a (z,p z ;c±, p) is strongly peaked around the char- 
acteristic momentum p z = p z and therefore Sj, should be 
a slowly varying function of ^. This will form the basis of 
our approximations below. For the time being we keep 
the equations exact. 

The dynamics for Sb follows from (|101f> . 



-s b (2,C; q,p) = 



-s b (z,c, q,p) 



— 5 6 (z,C;q, p) 
dz 



coh. 



(105) 



The coherent part reads 



—S b (z,C;q,p) 



coh. 



i — 7- / dps exp (i(ps/K) 



(106) 

S a (z,p z +Pd~;q,p) 



Pzh J Pz +PS 

= - 1 ^-r Sb ( z 'C;q, p) 

p z n 

- / d W ex P (KPs/K) D ( rr ) S a (z,p z +ps;q,p) 

Pzn J \PzJ 

where we used p z j (p z + ps) = 1 + ~D {ps /Pz) with 

00 

n—1 

Formally, Equation 1)103(1 allows to introduce a differen- 
tial operator for Sb, 



-s b (z,C;q, p) 



.q_p 



l+V 



coh. 

-m d 

Pz dC 



Sb (z,C; q,p) 



(108) 
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Since Sb was constructed to have a weak dependence on 
£ we expect that the expansion l(107|) can be relied on 
at least in an asymptotic sense. For the second term in 
((11)5(1 one obtains in a similar way 



— 5 b (z,C;q,p) 



incoh. 



m f 

-=*■ / d P ' z dp' exp(Kp' z /h)j(p'+p' z e z ) 

Pz J 

5 6 (z,C;q,p-p')- (109) 



ih d 

Pz dC 



This integro-differential equation can be further simpli- 
fied by separating off the solution of the coherent part 
(|108fl for a vanishing (^-dependence and by a Fourier 
transformation that removes the convolution in (|109fl . 
This is done by the introduction of a third and final aux- 
iliary function, 



Sc (z, C; q, p) := J dp exp (i^ 



q p vp 

z — 1 — - — 



x 5 6 (z,C;q,p) . (HO) 
It reads in terms of the Wigner function 
1 



s c (z, C; q, p) = 



drdp dp z 



x exp [ — i 
P 



(2nh) 
rq+pp . q p 

h p z h 



z + i( 



Pz - Pz 



X 



^- [ dt W(r + ze z ,p+p z e z ;t). (Ill) 
m p J 



If one knows the function S c the interference pattern is 
immediately obtained since 



and heading to the right, 

lim W(r + ze z , p + p z e z ; t) — for z < 0. (115) 

t — >oo 

With these conditions the initial function is obtained 
from i|lll|) by assuming that the Wigner function evolves 
freely (without decoherence) until it reaches the bound- 
ary to the decoherence region, z — 0. It reads 

S* c (0, C; q, p) = ——2 
(27m) 

[a a a /'.r q + P P, -,Pz-Pz 
x / drdp dp z exp — i h iQ — 



xf dz W (r + ze z + —p,p + p z e z ;o) . (116) 
J-oo V Pz J 

After solving the differential equation 1113|1 the resulting, 
possibly blurred, interference pattern Q z (r) is obtained 
by taking the inverse Fourier transform of 1112(1 



(117) 



Q (z; r) = J dq^ h S c (z, 0; q, j-q^j . 

The evolution equation ((113(1 shows clearly the hierarchy 
of decoherence terms involved in the dynamics. If one ne- 
glects the right hand side of ((113(1 altogether one obtains 
the diffraction pattern in the paraxial approximation, 



QlT a 0) (z;q) = S c [0,0;q, 



Pz 



(118) 



The first term in (|113fl describes the effect of decoherence 
in the paraxial approximation, the second term gives the 
corrections to the propagation beyond the paraxial ap- 
proximation, and the third term describes the modifica- 
tion of the decoherence due to those corrections. 



Q{z;q) = S c ( z,0; q, — q 



(112) 



as follows from 1104(1 . The evolution equation of S c is 
obtained from ((105(1 . It is now a differential equation, 



— s c (z,(; q,p) 



-^-i — z - p + Qe z ) S c 0, C; q, p) 

Pz \Pz 



h d 



T-T>[ — — ) V p S c (z,(;q,p) 



+ - q 

P 

-^l(^z- P + Ce z )v(^)Sc{zX^p). 
Pz \Pz J \Pz dQJ 

(H3) 

In order to find the initial function S c (0, £; q> p) we take 
into account that the initial state is localized in the left 
half space, 

W{r + ze z ,p+p z e z ;Q) = for z > 0, (114) 



B. Decoherence in the paraxial approximation 

In the simplest approximation we neglect the correc- 
tions in i|113|) due to the 2?-terms. In this case (11 1311 can 
be immediately integrated, 

^/V^-P + C^V 

Pz Jo \Pz J 



s c (z,C;q>p) = ex P 

x5 c (o,C;q,p) 



(119) 



It follows from 1(112(1 that the resulting interference pat- 
tern is characterized by 



' para 



(z;q) =exp 






. Pz Jq 


x Q 


(7=0) / \ 
)ara V^; Hy 

[-M 


= exp 


xQ(7 a =°) (z;q) 



(z' - z)q 



dz' 



'(t'-t)q 



dt' 



(120) 
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with t := zm p /p z . Using l|118|) and (|99|l we obtain the 
final pattern corresponding to a solution of the master 



equation l|93|l in paraxial approximation. 



Q<fit 0) M = J dr'dpdq^-^ exp (i 

x [ dp z [ dz'W (r' + — p + z'e z ,p+p z e z ;0 

J J-oo \ Pz 



Pz Jo V Pz 



(121) 



r 



With this result it is easy to see that the stationary 
treatment of decoherence in Sect. HVl is equivalent to the 
dynamic approach in the present section. To facilitate 
the comparison we treat the present problem with the 
method of Sect. IIVI Take the beam to be in a nontrivial 
stationary state at z = 0, 

Wbcam (r, p + p z e z ) = g (p z ) w (0; r, p) . (122) 

In the case of coherent evolution (JSJ) the interference pat- 
tern reads then 

Q(z;r)cx / dq e^ h [ dp z {jp z )w (z,p z ;q) 



with 

w(z,p z ;q) 



(123) 



1 



— p ) q/h 



2nny • dpdr exp 1 _i 1 1 

xw(0;r,p). (124) 

Using the same procedure as in Sect. IIVI one finds how 
decoherence events that take place at a constant rate 
R(z) in (z;z + dz) will modify the pattern l|123|) . The 
result is given by the expression in (|123|l if the w (z; q) 
are replaced by 



w(z,p z ;q) =cxp - / R(z') 



1 — rj I q I dz' 



xw{z,p z ;q), (125) 

where 77 is the corresponding decoherence function. 
Clearly, l|125|) is the analog of l|12(J|) in the case of a sta- 
tionary description. The only difference is the appear- 
ance in (|120fl of p z instead of p z , which occurs because 
the additional assumption of a strongly peaked longitu- 
dinal velocity distribution was necessary in the time de- 
pendent calculation. The strong similarity between the 
results (|125|) and l|120|) shows that the treatment of de- 
coherence in Sect. IIVI is indeed equivalent to solving the 
master equation in paraxial approximation. 



VII. CONCLUSIONS 

In this article we presented an analysis of Talbot-Lau 
matter wave interference that provides a quantitative 



prediction of the effects encountered in the experimental 
realization. It was shown that by describing the station- 
ary beam in terms of the Wigner function both the inter- 
action with the grating forces and the effects of marko- 
vian decoherence can be incorporated analytically. In 
addition, the formulation allows one to distinguish un- 
ambiguously the quantum phenomena from the effects of 
classical mechanics. 

Recently, our theory was successfully app lied to de- 
scribe experiments with large molecules |l2l Il3l IbH ITU 
l55l ]. Correspondingly, the discussion of decoherence ef- 
fects in the present article was confined to those mecha- 
nisms most relevant in the interference of fullerenes and 
biomolecules. Indeed, the interaction with gas molecules 
and the emission of heat radiation are expected to be rel- 
evant sources of decoherence for all large particles. Our 
formulation applies immediately to those since the bulk 
properties of the particles, such as the polarizability or 
the absorption cross section, were used to describe the 
environmental coupling. 

Other decoherence effects might become relevant as the 
particles increase further in complexity. In particular, 
those couplings that entangle the center-of mass motion 
with the rotation of the particle or with its internal de- 
grees of freedom become sources of decoherence. In these 
cases the observed loss of visibility, which is inevitable if 
the detection is insensitive to the relative coordinates, 
can be calculated with the same approach as discussed 
above. 

The grating interaction will also require a more re- 
fined treatment at some point. The eikonal approxima- 
tion ceases to be valid for particles of increasing size, be- 
cause they interact stronger and at the same time they 
will have a longer interaction time. A more careful eval- 
uation of the propagation through the grating will be 
needed in those cases. 

A final remark concerns the ease of incorporating de- 
coherence effects in the present formulation of matter 
wave interference. It draws heavily on the fact that one 
is able to separate the rate of decoherence events from 
their effect. It seems that this approach, that avoids the 
solution of a master equation in time, can be a transpar- 
ent way of treating markovian dynamics. It is vindicated 
by the comparison with the more conventional solution 
of a corresponding markovian master equation that was 
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presented in the last section of this article. 
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